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Abstract 



A sharp TH1/TH2 dichotomy has often been used to define the 
effects of cytokines on autoimmune diseases. However contradictory 
results in recent research indicate that the situation may be more com- 
plex. We build here a simple mathematical model aimed at settling 
the contradictions. The model is based on a neural network paradigm, 
and is applied using Partial Differential Equations (PDE). We show 
here that a TH1/TH2 paradigm is only an external view of a complex 
multivariate system. 
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1 Introduction 



Immunologists generally agree that in the course of many autoimmune dis- 
eases, TH1 type cells have a damaging effect, whereas TH2 type cells have 
a beneficial effect ]TJ. Thus the disease is aggravated by a high ratio of 
THf/TH2 T cells and cytokines. This view is contradicted by many exper- 
imental results showing an inverse effect |2[]: TH2 type cells may enhance 
disease, if administrated at a certain time ||, and TH1 type cells can stop 
the evolution of the disease, if given at the appropriate time [|, |J. Other 
works show that gene knock-out mice with no INF 7 or TNF a still man- 
age to produce the pathological features of a THl-dependent autoimmune 
disease |6|, [<]]. Based on results from EAE (experimental autoimmune en- 
cephalomyelitis) H and insulin dependent diabetes melitus (IDDM) |J [TIJ 
in Non Obese Diabetic (NOD) mice ||11|| , we propose a larger model explain- 
ing both the general agreement and the contradictory results. We suggest 
that the TH1 and TH2 cells are not the cause of the disease, but a marker 
for a more general steady state. To settle the contradictions, we use a nu- 



merical model and general mechanisms of safe information transmission (12 
to explain the discrepancy in the case of gene knock-out mice. We show the 
explanatory power of our model and its ability to produce the known features 
of TH1 type diseases. 



2 Model 



Our model is based on a neural network paradigm fL3 |, where the overall 
situation is not described by the properties of any one cell type, but by 
the steady state of the network. We use PDE in order to get a continuous 
description of the network. Each cell in the network represents a class of cells 
with the same general properties, at least from the most relevant aspect of the 
THI/TH2 dichotomy. The model is an ad hoc description of the dynamics, 
but most of the links we use have some experimental backing. We built the 
model with the smallest number of cell types needed to explain qualitatively 
the experimental results. 
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2.1 Experimental results 

EAE is an autoimmune disease studied mainly in rats and in mice, and is 



used as a model for human multiple sclerosis |14| [L5| . In EAE, myelin in the 
CNS becomes inflamed inducing paralysis and even death fTE] . This disease 
is classically defined by a TH1/TH2 dichotomy in which TH1 T cells induce 
the disease whereas TH2 cells prevent it [|TJ. 
Classical EAE experiments show that: 

• EAE can be induced either by adoptive transfer of CD4 T cells from a 
sick rodent, or by active immunization by an injection of a myelin anti- 
gen in adjuvant to a healthy rodent [ I7| . We will show that these two 



procedures lead, in fact, to two separate pathways of disease induction. 



Administration of free antigen arrests EAE [|TBJ. We propose to show 
that antigen-induced inhibition could be due to the effect of regulatory 
anti-idiotypic (anti-id) suppressor cells. 

The disease can be prevented by the continuous administration of IL-4, 
which is a TH2 type cytokine thought to push the immune system into 
a healthy TH2 dominated state REF FROM EAE|T§. But IL-4 can 
also aggravate disease (see below). 

However, recent research results appear to contradict the classical TH1 /TH2 
hypothesis: 

• A small amount of IL-4 administrated during disease induction can 
block the suppressive effect of free antigen. p(J This paradox will be 



explained by the effect of TH2 type cytokines on anti-id regulatory 
cells. 

• A short administration of a TH1 cytokine can heal the disease 0, |5|. 
This effect too can be attributed to an anti-id regulatory cell. 

• Gene knock-out mice, lacking one or some of the TH1 cytokines, can 
still develop EAE or IDDM || [7j. This paradox will be explained by 
the nature of regulation. 

These facts regarding EAE are combined in our analysis with results 
regarding the structure of id - anti-id cell interactions and the cytokine profile 
expressed during the evolution of IDDM in NOD mice |], [1(| , and in the 



mouse model of EAE || 21]. 
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2.2 Cells and Cytokine Types 



Our model is composed of five cell types, three cytokine types and the antigen 
(Table 1). The first cell type is the naive CD4 T cell that can differentiate 



into the second and third cell types: TH1 and TH2 type CD4 T cells |22 
The anti-id cells are grouped as a single cell type, the fourth cell type. A more 
detailed model should use naive anti-id cells and diverse anti-id cells secreting 
THf type cytokines and TH2 type cytokines. For the sake of simplicity, 
however, we group all anti-id cells into one cell type, and assume they have 
the same effect. Since we know from IDDM and EAE that anti-id cells secrete 
more TH1 cytokines than TH2 [TO, 0] , we designate the anti-id population 



as TH1 type cells. The CD4 T id cells can present themselves directly to 



the anti-id cells [10 . The fifth type of cells are macrophages, which play 



additional roles as APC and as cytokine-secreting cells. The antigen 
represents the quantity of antigen presented to the CD4 T cells. We assume 
that there is a large enough number of APCs so that the CD4 T cells can 
see the antigen. 

The three types of cytokines are: 

• TH1 type cytokines that enhance the differentiation of naive cells into 
TH1 type cells. 

• TH2 type cytokines that enhance the development of TH2 type cells. 

• Cytokines secreted by macrophages that enhance the proliferation and 
differentiation of TH1 cells [p4j] . We shall ignore the effect of TH2 type 



cytokines secreted by macrophages (like IL-10) and the effects of B 
cells. We will group together all the macrophage-derived cytokines as 
type C cytokines. 



2.3 Model equations. 

We compute the time derivatives of each cell type concentration Xi and of 
the cytokine concentration yi, assuming both are affected by all other cells 
and cytokines. 



Xj — ^iWijXiXj -\- YHtXX.^YjW\ijX or ig.i n iXj , 0) diXi (1) 

y\ = max(T l w 2 i j x j ,0) - (2) 
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Idiotypic TH1 cells 


A 


3 


Idiotypic TH2 cells 


B 


4 


Anti-id cells 


anti-id 


5 


Macrophages 


M$ 


6 


TH1 cytokines 


Cyt A 


7 


TH2 cytokines 


Cyt B 


8 


Cytokines secreted by macrophages 


Cyt C 



Table 1: Cell types and cytokines. The cytokines are labeled according to 
the cell type producing them. We grouped cells as follows: TH1 id cells can 
include all types of lymphocytes producing TH1 cytokines. Anti-id cells are 
all cell types specifically reactive to the id cells. They can either be cytotoxic 
cells, or limit the proliferation rate of the id cells. 



Xorigini is the origin of the cell number i (for example, the TH1 cells 
originate in the naive cells, the naive cells originate in themselves). The 
other cells and cytokines can either induce or inhibit the proliferation of this 
(xi) cell type. The cytokines originate from the cells that secrete them, so 
that the right side of eq' 2 is not multiplied by the cytokine concentration. 
The secretion rate of the cytokines (yi) by cells (xj) is linear to the other 
cell and cytokine concentrations. For example, if TH2 type cytokines inhibit 
the production of TH1 cytokines, then we use a term in the equation for the 
TH1 cytokines —W2ijx(th,2), and we do not multiply the term by the TH1 
cytokine concentration. 

A cell's proliferation rate cannot be negative. Thus we define the prolif- 
eration rate of any cell to be greater than 0. But of course, the natural death 
rate (—diXi) can lead to a decrease in the cell concentration. T he cytokine 
production rate (y'j) is assumed not to depend on the cytokine concentration, 
but on the number of cells secreting the cytokine. The equations are solved 
numerically using a second order precision scheme |^5| . 

We start every computation at a steady state 



x 



i = 0,<4 = (3) 



Since the time constants in this process are short, we could assume that 
the system had already reached a steady state before the appearance of the 
autoimmune disease. 
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Macrophages do not reproduce during the response to the antigen but, 
rather, are attracted to the affected region [p6j] , so that, like the cytokines,their 



growth rate is not multiplied by their own concentration. 

The differences between possible models are in the assigned weights Wij, wuj 
and W2ij in eq' 1 and 2. But we have tried to use the most uniform definition 
of weights, and all the weights are of the same order of magnitude. If we had 
no special reason to believe that two weights of the same type were different, 
we fixed them at the same value (the effect of TH1 cells on TH1 cytokines 
and the effect of TH2 cells on TH2 cytokines, for example). This was done 
in order to limit the size of the parameter space. 

2.4 Model description 

The structure of the model is simply a mathematical description of the find- 
ings in IDDM and EAE. The basic part of the model is the opposition between 
id cells and anti-id cells. The CD4 are the id cells and the naive CD4 popula- 
tion proliferates proportionally to the antigen concentration. The suppressors 
are the anti-id cells, and they proliferate proportionally to the concentration 
of id cells (Figure 1A). 

The second division of this model is into TH1 and TH2 id cells. Both 
the TH1 and TH2 cells originate from the naive id cell. The TH1 cell type is 
associated with a TH1 cytokine type and the TH2 cell type is associated with 
a TH2 cytokine type. We assume that each cell enhances its own cytokine 
type and inhibits the production of the opposing type. Each cytokine 



encourages the production of the appropriate cell type and inhibits the other 
type (Figure IB). 

The macrophages have a double effect: They serve as APCs and activate 
naive CD4 cells |23|, and they secrete cytokines that enhance T cell differen- 
tiation into TH1 type cells p8f (and inhibit differentiation into TH2 cells). 
The macrophages are influenced by the TH1 cytokines, which activate them 
to migrate into the inflammatory site. Since macrophages respond to TH1 
cytokine and are inhibited by TH2 cytokines, changes in the macrophage 
concentration depend on the differences between the concentrations of TH1 
and TH2 cytokine types. 

The anti-id cells are of a TH1 type. This assumption is based on data 
from EAE and IDDM |I| ||. We thus postulate that the anti-id cells are 
inhibited by TH2 type cytokines. And finally, we assume that there is a basic 
concentration and differentiation rate of CD4 cells into TH1 and TH2; but 
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Figure 1: The basic features of the model. 1 A- Regulation of the id cell 
population by anti-id cells. IB - Competition between TH1 T cells and TH2 
T cells. Full lines designate activation. Broken lines designate suppression. 



the natural situation (when there are no cytokines involved) tends toward the 
TH2 dominated state (a healthy situation). The above model is described in 
Figure 2 



2.5 Model mechanism 

We have two basic steady states (Figure 3), and all the cells and cytokines 
administrated only move the model from one steady state to the other, and 
back. The basic situation is a TH2 type steady state, which is assumed 
to be healthy, whereas the TH1 steady state leads to autoimmune disease. 
In contrast to the general view, the TH2 steady state represents a global 
situation of the system, and not the concentration of a single cell or cytokine. 
Note that the absence of disease does not mean that the system is at rest; the 
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Figure 2: Flow chart of the model. The model is based on a TH1/TH2 
balance. If the TH1 concentration is too high the target organ is attacked. If 
there is a high TH2 concentration, the target organ remains healthy. Dotted 
lines indicate inhibition; full lines indicate stimulation. There is also an 
inhibition of TH2 on the production of TH1 and vice versa, but these lines 
were not drawn in order not to overload the graph. 
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TH2 steady state is associated with relatively low concentrations of reacting 
cells. Nevertheless the id and ant-id cells are still reacting at the same rate, 
but at a lower concentration. In other words the healthy state merely looks 
like it is at rest; on the contrary health like disease is a state of interactions 
[29 . 

When we administer cells or antigen and adjuvant, we push the system 
toward the TH1 state (Figure 4). Following administration of T cells, which 
are primarily TH1 type cells, the number of TH1 cells is high enough to 
start the production of TH1 cytokines that inhibit the TH2 cytokines. The 
system then falls into the TH1 steady state, which is the diseased state. 
Another way to induce disease is to immunize the animal with antigen and 
adjuvant. We propose that the adjuvant stimulates macrophages to produce 
C type cytokines. The antigen together with the cytokines induces CD4 T 
cell proliferation and differentiation into TH1 cells, bringing the system again 
to a TH1 type steady state. 

One way to heal the disease is to administer free antigen (Figure 5) in 
an amount that can raise the total number of naive CD4 cells that become 
TH1 CD4 cells. The resulting high concentration of CD4 cells will activate 
an increase in the number of anti-id cells. At this point, we assume that the 
lifespan of the anti-id cells is longer than that of the CD4 cells; the number of 
anti-id cells will remain high even after the number of CD4 id cells decreases 
(due to the fall in antigenic stimulation). Once the total number of CD4 cells 
is below a certain threshold, the system returns by itself to the healthy TH2 
situation. The basic situation is the TH2 situation. 

The return to a TH2 state, following antigen administration, can be para- 
doxically inhibited by the administration of IL-4 with the antigen (Figure 6). 
The reason is as follows: When we administer IL-4 (which is represented 
by type B cytokines), we inhibit the production of anti-id cells. In the ab- 
sence of a sufficient number of anti-id cells, the number of CD4 effectors cells 
remains above the threshold and the system stays in a TH1 state. If, how- 
ever, the quantity of IL-4 is high enough to overcome TH1 dominance and 
starts activating TH2 T cells directly, then the TH1 cell concentration will 
decrease due to the production of TH2 cytokines by the newly formed TH2 
cells, despite the low number of anti-id cells. 

The administration of anti-IL-4 antibody with the antigen can paradoxi- 
cally act in a manner similar to the administration of IL-4 with the antigen. 
Note that in both cases the inhibition of disease is weaker than the effect of 
administrating free antigen. The administration of anti-IL-4 antibody has 
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TH1 vs TH2 concentration as a function of time. 
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3B TH1 vs TH2 concentration as a function of time. 
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Figure 3: The basic steady states. The first steady state is TH1 dominated 
(3A) and the second is TH2 dominated (3J3). The model is such that, starting 
from a naive state, the system equilibrates into a TH2 steady state. 3A and 
3B show concentrations in arbitrary units. Both 3A and 3B are composed of 
4 graphs. The evolution in time of the TH1 and TH2 T cells concentration 
is depicted in the upper left drawing^ The upper right drawing shows the 
cytokines concentration as a function of time. The drawing in the second 
row to the left shows the average concentration of id cells (TH1, TH2 and 
naive cells), the concentration of anti-id cells and the concentration of free 
antigen. The drawing at the right of the second row show the evolution in 
time of the id and the anti-id cell concentration. Note that the TH2 steady 
state (3B) manifests a lower concentration of both id and anti-id cells. 
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TH1 vs TH2 concentration as a function ot time. 
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Figure 4: The Transition from TH2 (healthy) steady state to a TH1 (disease) 
steady state. 4A shows the effect of adminstrating (arrow) TH1 cells. 4B 
shows the effect of immunization with antigen and adjuvant (arrow). Note 
that the TH1 steady state is associated with enhancement of the anti-id cells. 
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Figure 5: Healing with free antigen. Administrating free antigen raises the 
TH1 cell concentration, and then stimulates anti-id cells. Folowing it the 
anti-id cells concentration rises. The anti-id in turn supress the system and 
it returns back to its natural TH2 steady state. The notation used in this 
figure is identical to the notation used in figure 3. 
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TH1 vs TH2 concentration as a function of time. 
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Figure 6: Inhibition of the healing effect of free antigen by administration 
of IL4. IL4 has a double role in this case. At the first stage it reduces the 
number of TH1 id cells, and at the second stage it reduces the number of 
anti-id cells, so that their effect is weakened. The notation used in this figure 
is identical to the notation used in figure 3. 
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no effect on the system at first, since the concentration of IL-4 is low. The 
effect of the antigen is then as described at the beginning of the section; the 
antigen induces a sharp rise in the concentration of id CD4 cells, and follow- 
ing the id CD4 cells, the anti-id cell concentration raises. These anti-id cells 
produce a long term decrease in the concentration of the CD4 id cells. At 
a later stage, when the concentration of CD4 id cells has decreased to a low 
level, the effect of the anti-IL-4 antibody becomes important. The level of 
id CD4 TH1 type cells at this stage is low, and the system should now start 
to produce TH2 type cells and cytokines. However the anti-IL-4 antibody 
now inhibits the production of IL-4, and the system has a lower probability 
to return to a TH2 state. 

Just as the administration of TH2 cytokines can inhibit the transition 
back to a TH2 state, the administration of TH1 cytokines can rise the prob- 
ability of such a transition. The continuous administration of TH1 cytokines 
will of course keep the cytokines in the TH1 profile and induce or enhance 
the disease. Administration of the TH1 cytokines at the onset of the disease 
(Figure 7) will lead to an increase in the production rate of TH1 id CD4 
cells; the high concentration of TH1 cells will raise the concentration of the 
anti-id cells, which will later feed back to decrease the concentration of CD4 
id cells. If the CD4 id cell concentration decreases to a low level, the sys- 
tem will return to its natural steady state, which is a TH2 steady state, as 
described above. 

The administration of some cytokines has been shown to have no observ- 
able effect on the system. This does not mean that these cytokines have no 
important role in the dynamics of this system. In fact the lack of an observ- 
able effect can be due to a high endogenous concentration of the cytokines; 
adding additional amounts cannot enhance the effect. 

Administration of a cytokine may have no effect if the cytokine concen- 
tration is close to, or above its threshold for causing an effect. For example, 
administration of IL-12 at the onset of disease may have no effect if the IL-12 
effect is already at its maximum. Paradoxically in this case, the administra- 
tion of a low dose of anti-IL-12 antibody will also have no effect, since the 
precise concentration of the molecule is not important. A large change in the 
concentration of the molecule leads to a small change in its effects. This can 
be understood if the effect is bell shaped as a function of the concentration. 
At low concentrations, a small change in the concentration leads to a large 
change in the function value. However at the top of the bell, a large change 
in the concentration leads to very small change in the function value. 
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TH1 vs TH2 concentration as a function of time. 
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Figure 7: Healing by administrating TH1 cytokines. The TH1 Cytokines 
raise the TH1 cell concentration, inducing a rise in the concentration of anti- 
id cells. The anti-id cells lower the concentration of the id cells, returning the 
system back to its normal (healthy) TH2 steady state. The notation used in 
this figure is identical to the notation used in figure 3. 
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2.6 Redundancy and gene knock-out mice 



Gene knock-out mice lacking some of the TH1 cytokines can still develop 
a THl-dependent disease. However this does not contradict the model we 
presented, which describes the evolution of the disease as the transition to 
a TH1 type cytokine profile. The cytokines in this model are described 
as a group with the same approximate action. Even if such a group lacks 
one of its components, the effect of the group as a whole is not altered. 
The TH1/TH2 phenotype does not depend on a single cytokine, since the 
concentration of each one of the TH1 and TH2 cytokines can change by a 
large number of mechanisms. The immune system cannot afford to change 
its general phenotype each time an external factor modifies the concentration 
of a single cytokine. The immune system achieves its robustness to external 
changes in the cytokine profile through redundancy ||29|| . A large number of 
cytokines can each have common effects. Thus a certain result can occur 
even if not all the cytokines leading to the effect are expressed. 

Indeed, the fact that the body uses many redundant cytokines can have a 
larger meaning. Consider that the effects mediated by cytokines have at least 
two functions. The first function is to maintain a healthy steady state, and 
the second function relates to the fine details of the steady state. In normal 
behavior, the changes are felt at the level of the details, but in extreme 
situations, like the response to an external pathogen, the body will try to 
change its total steady state and produce a fast and strong response. 



3 Discussion 

3.1 Regulation of autoimmunity 

This paper presents a simplified model of the kinetics of the concentrations of 
the main cell and cytokine types in the evolution of a TH1 type autoimmune 
disease. It shows that even simple dynamics can explain many of the complex 
features in the evolution of the response. We have grouped the major cells 
and cytokines into single types, although there are obviously differences in 
the roles taken by each one of the cells we have grouped together. This 
model ignores the effects of localization and of tissue barriers (for example, 
the blood/brain barrier in EAE) in inhibiting the passage of some of cells and 
allowing others. We did not take into account the effects of the tissue cells 
(for example astrocytes in EAE, or islet cells in IDDM), of the cytokines they 
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secrete, and of B cells and antibodies. Some major mechanisms were ignored 
including anergy and the release of free antigen due to cell destruction. 

All these simplifications render this a general feature model rather than 
a precise quantitative description of the cell concentrations. Our model is 
built on a double system of positive and negative feedback between id and 
anti-id cells, and between TH1 and TH2 states. Such a structure could 
be the basic feature of the selection between TH1 and TH2 for a group of 
important antigens in the body. This model seems to describe the situation 
both in EAE and in IDDM, and we have no reason to believe that it is 
not general. Negative feedback keeps the general concentrations of all kinds 
of cells within reasonable limits, and positive feedback enables the system 
to mount a fast enough response to an antigen when needed. With this 
structure, it is possible to maintain a population of T cells responding to 
self-antigens of the body, keeping them in a TH2 state where they may be 
restrained. 

3.2 Models 

This kind of model shows the complex dynamics attained even by a small 
system containing only 5 cell types and 3 cytokine types. It shows that 
even small systems need a kinetic description, and the system cannot be 
understood by an overall estimation of the effect of one cell type on the 
others. Such models must be robust to changes in parameters, since it is 
obvious that the interactions between the cells in the body are not precisely 
fixed and that the cells are constantly subjected to chemical noise impinging 
on them. Nevertheless, the body must keep its general steady state, although 
the precise quantity of any cell is not very important. 

3.3 Feedback loops 

The present type of modeling shifts our attention from the numbers of each 
cell type to the global state of the immune system, which most of the time 
should be in some quasi-equilibrium between many opposing influences. The 
way to analyze such a global state is by positive and negative feedback loops 
and the interactions between them. According to this formalism, our model 
can be described as a large negative feedback loop, built of 4 internal loops : 

• The id-anti-id cell negative feedback loop, in which are nested three 
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positive and one negative feedback loop. These include: 

• The TH1 and TH2 cells each of which feedsback positively on itself. 

• The TH1 and macrophage positive feedback. 

• The TH2 and anti-id cell negative feedback. 

In conclusion, the description of disease evolution can be based on a global 
kinetic view, and the subsystems are not cell types but feedback loops, (figure 
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Figure 8: A schematic description of the five feedback loops. 8A Shows the 
negative feedback loop between id an^othe anti-id cells. 8B and 8C show the 
positive feedback of the TH1 and TH2 cell and cytokines on themselves. 8D 
shows the positive feedback of the macrophages and the TH1 id cells and 
macrophages. 8E shows the negative feedback between the TH2 cytokines 
and the anti-id cells. 
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